COUPLING IMPEDANCES OF SMALL DISCONTINUITIES: 
DEPENDENCE ON BEAM VELOCITY 
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The beam coupling impedances of small discontinuities of an accelerator vacuum chamber have 
been calculated [e.g., S.S. Kurennoy, R.L. Gluckstern, and G.V. Stupakov, Phys. Rev. E 52, 4354 
(1995)] for ultrarelativistic beams using the Bethe diffraction theory. Here we extend the results 
to an arbitrary beam velocity. The vacuum chamber is assumed to have an arbitrary, but uniform 
along the beam path, cross section. The longitudinal and transverse coupling impedances are 
derived in terms of series over cross-section eigenfunctions, while the discontinuity shape enters 
via its polarizabilities. Simple explicit formulas for two important particular cases — circular and 
rectangular chamber cross sections — are presented. The impedance dependence on the beam 
velocity exhibits some unusual features: for example, the reactive impedance, which dominates in 
the ultrarelativistic limit, can vanish at a certain beam velocity, or its magnitude can exceed the 
ultrarelativistic value many times. In addition, we demonstrate that the same technique, the field 
expansion into a series of cross-section eigenfunctions, is convenient for calculating the space-charge 
impedance of uniform beam pipes with arbitrary cross section. 

PACS numbers: 41.75.-i,41.20.-q 



I. INTRODUCTION 

A common tendency in design of modern accelerators is 
to minimize beam-chamber coupling impedances to avoid 
beam instabilities and reduce heating. Even contribu- 
tions from tiny discontinuities like pumping holes have 
to be accounted for because of their typically large num- 
bers. Direct numerical methods encounter difficulties 
in calculating the impedances of small discontinuities. 
The difficulties of usual time-domain methods — com- 
puting wake potentials and then finding the impedances 
as their Fourier transforms — are mostly technical for 
ultrarelativistic beams: one has to apply very fine steps 
in wake computations to resolve small obstacles. How- 
ever, for non-ultrarelativistic beams the usual numerical 
approach fails due to difficulty in implementing proper 
boundary conditions at the open ends of the beam pipe. 
This situation makes analytical methods especially im- 
portant. A general analytical approach for calculating 
the beam coupling impedances of small discontinuities on 
the walls of an accelerator vacuum chamber has been de- 
veloped in [10 for ultrarelativistic beams. The method 
is based on the Bethe theory of diffraction by small holes 
, according to which the fields diffracted by a hole can 
be found as those radiated by effective electric and mag- 
netic dipoles. 

The Bethe idea of effective dipoles developed in 1944 
was first used to calculate the coupling impedances of 
pumping holes in a circular waveguide for the case of an 
ultrarelativistic beam in 1992 ^ Bethe's theory is 
applicable for wavelengths large compared to the typical 
hole size h, and it can be used for impedance calcula- 
tions when h is much smaller than the typical dimension 
b of the chamber cross section. Within this limitation, 
the theory was applied to vacuum chambers with an ar- 
bitrary simply-connected cross section in ,2j], again for 
relativistic beams. The imaginary part of the impedance 



was shown to be proportional to the difference of hole po- 
larizabilities (ip — x), where the magnetic susceptibility ip 
and the electric polarizability x are both small compared 
to b 3 . From considerations of the energy radiated into the 
chamber and through the hole, the real part of the hole 
impedance comes out to be proportional to (ip 2 + \ 2 ), 
being usually much smaller than the reactance. These 
results are not restricted to small holes, but remain valid 
for other small discontinuities like posts, masks, or irises, 
because the idea of effective polarizabilities works equally 
well in those cases, as was demonstrated in 0,0- A more 
refined theory that takes into account the reaction of ra- 
diated waves back on the hole was developed in In 
the latter approach, the beam coupling impedances of a 
discontinuity come out as perturbative series in polar- 
izabilities, more exactly, in small parameters ip/b 3 and 
x/b 3 : the reactive impedance is the first order effect, 
while the real part of the impedance has the second or- 
der. Moreover, the theory [l| contains non-perturbative 
effects: it gives the trapped modes near the pipe cut-off 
frequencies due to some small discontinuities. This effect 
was discovered earlier using a different method 

In the present article we extend the analytical ap- 
proach P, 13 to the case of non-ultrarelativistic beams. 
The beam coupling impedances of a small discontinu- 
ity on the walls of a vacuum chamber with any simply- 
connected cross section are derived for an arbitrary beam 
velocity v — (3c. The theory gives analytical expressions 
for the longitudinal and transverse coupling impedances 
in terms of series over cross-section eigenfunctions. The 
shape of the discontinuity enters via its electric and mag- 
netic polarizabilities. Previous results for ultrarelativis- 
tic beams are naturally reproduced in the limit of f3 — > 1 . 
Simple explicit expressions for the impedances of non- 
ultrarelativistic beams are derived for circular and rect- 
angular chamber cross sections. 

One should mention a few earlier results related to 
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the subject. The coupling impedances of a small hole 
in a circular beam pipe were calculated @ for (3 < 1 
in the form of series involving roots of Bessel functions. 
Closed-form expressions of the hole impedances for non- 
relativistic beams in a cylindrical vacuum chamber with 
a circular cross section were obtained by Gluckstern et 
al [ToL ITU ] . The longitudinal impedance of a small round 
hole in a circular waveguide was derived in a closed form 
for a beam with (3 < 1 and a finite transverse size in • 
In a somewhat related paper ^3|> the bunch loss factors 
were investigated for non-ultrarelativistic beams. All the 
results above, however, have been restricted to the par- 
ticular case of an axisymmetric vacuum chamber. Our 
study treats a more general case of an arbitrary chamber 
cross section. 

In addition, we show that the same technique - 
based on field expansions in cross-section eigenfunctions 
- works very efficiently for deriving the space-charge 
impedance of a uniform waveguide having an arbitrary 
cross section. A general expression is obtained that in- 
cludes frequency corrections to and beam- velocity depen- 
dence of the space-charge ^-factor. In the particular case 
of a circular cross section, the results coincide with the 
known ones, e.g. For a rectangular chamber, the 

space-charge impedance is expressed in the form of a sim- 
ple series convenient for computation. 

The paper is organized as follows. A general anal- 
ysis of the beam fields and fields scattered by a small 
discontinuity in the vacuum chamber is given in Sec. II. 
Section III presents the derivation and results for the cou- 
pling impedances of small discontinuities. In Sec. IV the 
space-charge impedance of a homogeneous chamber with 
an arbitrary cross section is calculated. The eigenfunc- 
tions and some derivations for the two particular cases of 
the vacuum chamber — with a circular and rectangular 
cross section — are presented in the Appendices. 



II. FIELDS 

Let us consider an infinite cylindrical pipe with an ar- 
bitrary cross section 5" and perfectly conducting walls. 
The z axis is directed along the pipe axis, a small dis- 
continuity (e.g., a hole) is located in the cross section 
z = at the point (b, 0), and a typical hole size h satis- 
fies h <C b. The discontinuity is considered small when its 
size is much smaller than the wavelength of interest. To 
evaluate the coupling impedance one has to calculate the 
fields induced in the chamber by a given current. Con- 
sider a charge q that moves on or parallel to the cham- 
ber axis with velocity v — f3c. A rigid bunch is usually 
assumed for wake or impedance calculations. It is con- 
venient to choose a "pancake" charge distribution such 
that the charge and current densities are given by 

p(r,z;t) = qf(f)S(z~(3ct), 

j(r,z;t) = q[3cf{r)5{z- I3ct)z , (1) 



where f(r) is a normalized transverse charge distribu- 
tion, J s drf(r) = 1, S(x) is Dirac's delta function, and z 
is a unit vector in the z-direction. The distribution i[TJl 
includes two particular cases that are especially conve- 
nient in calculations: a point charge with the transverse 
offset s from the axis, when f(r) = 5(r — s), and a thin 
axisymmetric ring of radius a with f(f) — 8{r—a)/(2Tta). 
In frequency domain, omitting factor exp(— iut), the cur- 
rent density harmonic becomes 

j(r,z;u J ) = qf(r)e i » z /^z . (2) 

This harmonic corresponds to the space harmonic k = 
u>/ (f3c). From Maxwells's equations in frequency domain 
follow the wave equations for harmonics of the beam field 
components: 

(V 2 + d 2 ) E z + u 2 /c 2 E z = nup /(f3 2 7 2 )j z , 
(V 2 + d 2 )E ± +u 2 /c 2 E ± = l/{e (3c)V Jz , (3) 
(V 2 + d 2 ) H ± + u?lc 2 Ry_ = -(ixvjj,. 

Here all quantities are functions of (r, z;ui), V is the 
two-dimensional (2D) gradient in plane S, d z = d/dz, 
7=1/ y/l — ft 2 , and j z is defined from Eq. . The lon- 
gitudinal component of the magnetic field is not excited 
by the current 0J, so that H z = 0. 

Let us now introduce eigenvalues k 2 and orthonorma- 
lized eigenfunctions (EFs) e g (r) of the Dirichlet bound- 
ary problem in S: 

(V 2 + fc 2 )e s = 0; e 9 | as = 0, (4) 

where g — {n, m} is a generalized 2D index. Since EFs 
form a complete set in S, the transverse dependence of 
solutions to Eqs. © can be found as a series in EFs. Ob- 
viously, the solutions should depend on the longitudinal 
coordinate as exp(iujz/ (3c). Then the fields harmonics 
E, H produced by the charge distribution at the lo- 
cation (6, z) on the chamber wall without hole can be 
expressed in terms of EFs (@J as 

E u (b,z;uj) = Z H T (b,z;uj)/(3 (5) 




where it was convenient to introduce k = w/(/?7c). Here 
Zq = \J Mo/eo = 1207T Ohms is the impedance of free 
space, v means an outward normal unit vector, f is a 
unit vector tangent to the boundary dS of the chamber 
cross section S, V„ = V • v, and {u, f, z} form a right- 
handed basis. In Eq. © f g are the coefficients of EF 
expansion /(f) = J2 g fg e g( ff )'-> the Y are g iven b y fg = 
J s df f(f)e g (f). For the case of a point charge with the 
transverse offset s from the axis, we have f g = e g (s ). 
The eigenvalues and EFs for particular cross sections are 
given in the Appendix. 
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In a similar way, the longitudinal component of the 
electric field produced by the current Q as a function of 
the transverse coordinates f is 



E z (r, z;u>) = -i 



u Z q iuz/uw f g e g (r) 



13c /3 7 2 



^ k 2 



(6) 



Obviously, E z vanishes at the chamber wall because of 
the Dirichlet boundary conditions (BCs) for EFs (QJ . The 
transverse fields on the wall JSJ) are not zero since the 
series © includes the gradients of EFs and V ' v e g (b) ^ 0. 

Note that one could use EFs different from EFs |@J| 
to calculate the transverse magnetic field H±. A conve- 
nient choice would be to use EFs h g (f) of the Neumann 
boundary problem, which differs from J3J by the BCs, 
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0. However, we chose to relate H± to the 

transverse electric field Ej_ directly from the Maxwell 
equations. For the harmonic k = ui/(/3c), one gets 
Zqz x H±_ = —j3E±, from which the first line in Eqs. © 
follows immediately. 

One more remark about Eqs. JSJ is worthwhile. Since 
the lowest eigenvalue k g is of the order of 1/6, in the 
long-wavelength limit, nb = ujb/(f3jc) <C 1, the sums 
in Eqs. JH}-® become independent of frequency and of 
[3. The only remaining dependence on the beam veloc- 
ity E v oc 1/(3 seems counter-intuitive: one would expect 
H T oc [3. The apparent contradiction disappears as soon 
as we recall that Eqs. |J5J give the field harmonics with 
wavenumber k = w/(/3c), not the beam fields in time 
domain. 

Now we have to calculate the fields scattered into the 
vacuum chamber by the discontinuity. According to the 
Bethe theory, at distances / such that h -C I -C b, the 
fields radiated by the discontinuity (hole) into the pip e 
are equal to those produced by effective dipoles 
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P v = -xeo^/2; M T = ^ TT H^ +iP TZ H^)/2; 
M z = (V>»rfi?+V„fl£)/2, (7) 

where superscript '/i' means that the beam fields (JSJ are 
taken at the hole location (6,0). Polarizabilities are 
related to the effective ones a e ,a m used in 0, |TI| as 
a e = — x/2 and a m = ip/2, so that for a circular hole 
of radius h in a thin wall ip = 8h 3 /3 and \ = 4/i 3 /3 

0. In general, ip is a symmetric 2D-tensor, which can 
be diagonalized. If the discontinuity is symmetric, and 
its symmetry axis is parallel to z, the skew terms vanish, 

1. e. ip TZ = tp ZT = 0. In a more general case of a non-zero 
tilt angle a between the major symmetry axis and z, 



ijj TT = ip± cos 2 a + ipn sin 2 a , 

"0rz = ipzr = [ip\\ — ip±) sin a cos a , 

"0zz = ip± sin 2 a + ipn cos 2 a , 



(8) 



where ipu is the longitudinal magnetic susceptibility (for 
the external magnetic field along the major axis), and 
ip± is the transverse one (the field is transverse to the 



major axis of the hole). Formulas for polarizabilities of 
various discontinuities are collected in the handbook |l5j , 
see also Ref. [T^ . 

When the effective dipoles are obtained by substitut- 
ing beam fields (J5J into Eqs. Q, one can calculate the 
scattered fields as a sum of waveguide eigenmodes excited 
in the chamber by the dipoles, and find the impedance. 
This approach has been carried out for a circular pipe in 
4] , and for an arbitrary chamber in 0] , for ultrarelativis- 
tic beams. The fields radiated into the chamber by the 
effective dipoles Q can be found as a series in TM- and 
TE-eigenmodes [l4|: 



F = 



J2[AtF^ + 9(z) + A g F^-9(-z) 
g 



(9) 



where F means either E or H, g — {n, m} is a 2D index; 
superscripts '±' denote waves radiated respectively in the 
positive (+, z > 0) or negative (— , z < 0) direction, and 
9{z) is the Heaviside step function. The fields Fg of 
the gth. TM-eigenmode in Eq. © are expressed [Tj] in 
terms of EFs (QJ 



Ef = fc 2 e 3 exp(±r 3 z); Hj = ; 
Ef = ±r 9 Ve 9 exp(±r 9 z); 
Hf = iujeqz x Ve s exp(±r g z) , 



(10) 



where propagation factors T g — (fc 2 — oj 2 /c 2 ) 1 / 2 should be 
replaced by — i(3 g with [3 g = (lu 2 /c 2 — kg) 1 / 2 for uj/c > k g . 
For given values of dipoles Q the coefficients can be 
found 

Si 

using the Lorentz reciprocity theorem 



with 



Af = a g M T ± b g P v , 



.^L V e h - b = — 
2Y n kV vG ^ bg 2e k 2 



9 n g 



(11) 



(12) 



Here the EFs are taken at the hole location b on the wall, 
e h g =e g (b). 

Similarly, the fields F g ' of the TE g -eigenmode in 
Eq. |JSJ| are 



Hf = k' 2 h g exp(±r' g z); Ef=0; 

Hf = ±r;v/i s exp(±r^) ; 

Ef 1 — —iui[iQZ x Vh g exp(±T' g z) , 



(13) 



with propagation factors T' g = (k! 2 — lo 2 /c 2 ) 1 / 2 replaced 

by -i/3' g = -i(uj 2 /c 2 - k 12 ) 1 ' 2 when u/c > k' g . Here EFs 
h g satisfy the boundary problem (0J with the Neumann 



boundary condition, V^/i g |g5 

sponding eigenvalues, see in Appendix. The TE-mode 



0, and kg are corre- 
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excitation coefficients in the expansion © for the radi- 
ated fields are 



where 



Bf = ±c g M T + d g P v 



q g M z 



(14) 



(15) 



(b, 0), leads to 



Z(uj) = -iZ x 

c 2 

k 2 . + K 2 2-1 fc2 + K 2 



(18) 



It is convenient to introduce the following notation 



It \ — h"^ " eh g 

e v {f; K ) = -l^ w — 5 



(19) 



From Eqs. (f^Tf)- (|1 5|) . using the effective dipoles (J7J in- 
duced on the discontinuity by the beam fields (J5J , we find 
the fields scattered by the discontinuity into the vacuum 
chamber. Note that the dependence on the beam velocity 
enters explicitly only in Eq. (J5J. Now we proceed with 
the impedance calculation. 



III. BEAM COUPLING IMPEDANCE OF A 
SMALL DISCONTINUITY 

A. Longitudinal Impedance 

In a general case, the longitudinal impedance is defined 
(e.g., |l7l llSjl as an integral along the vacuum chamber 
of the (synchronous with the beam) harmonic of the lon- 
gitudinal electric field E z created in the chamber by the 
leading charge (QJ divided by the amplitude of the cor- 
responding harmonic of the beam current Jy. The field 
E z should be taken at the longitudinal position of a unit 
test charge that follows the leading one with the same ve- 
locity, and should also be integrated over the test-charge 
transverse distribution t(r). This leads to the formula 



Z(u) 



dze'^ / drt(f)E z (f,z;u), (16) 



where normalization J s dft(f) — 1 is assumed. Note 
that the longitudinal field E z includes both the beam 
field © and the scattered field © produced by the dis- 
continuity. The first one is related to the space-charge 
impedance and vanishes at 7 — > 00; we will consider it in 
Sec. IV. Here we include into E z only the discontinuity 
contribution. All z-dependence in Eq. J^J is in the expo- 
nents, so integrating over z is straightforward. For trans- 
verse integration, we expand the test-charge transverse 
distribution t(f) in EFs as t(r) = z\^ g t g s g (r), where 
t g = J s drt(r)e g (r). The result is 



Z(v) 



qc 



Z a M T 



IP* 

(3 £ 



^k 2 g 

9 9 



(17) 



where we again use variable k — to/ (Pjc) introduced ear- 
lier in Eq. ©. Substituting the dipole expressions (7J|, 
where the beam fields (jSJ) are taken at the hole location 



Comparing to Eqs. (JSJl, one can see that e v (f; n) is related 
to the transverse harmonics of the beam field at the hole 
location as 



E v (b,0;w)/3 _ H T {b,0;oj) 
Zoq q 



(20) 



As was already mentioned in Sect. II, the lowest eigen- 
value kg is of the order of 1/6. Therefore, for uib/(/3jc) = 
d « 1, the normalized transverse field Eq. 1)19(1 becomes 
frequency- and velocity- independent: 



e„(/;0) 



g 



fgkg V^Cg 



(21) 



The condition w6/(/?7c) <C 1 includes two cases: (i) ul- 
trarelativistic limit, 7 — > 00; and (ii) long- wavelength (or 
low-frequency) limit, when the wavelength A = 2tt /k — 
2tt[3c/ui is large compared to the typical cross-section size 
b. Equation l|21|) has a simple physical interpretation. It 
gives a solution for a 2D electrostatic field created on 
the chamber wall in the cross section S by a uniform in 
z charge, equal to e P er un it length of z, that has the 
transverse distribution f(r). From the Gauss law, e u (f) 
satisfies the normalization condition 



OS 



dl e„(/;0) = 1 



(22) 



where integration goes along the cross-section boundary 
dS. For a simple particular case of a circular cross section 
with radius b, and an axisymmetric charge distribution 
/(O = fi r ) — it includes an on-axis point charge, - 
Eq. (1221) by itself provides the solution,due to the prob- 
lem symmetry: e„(/;0) = 1/(2tt6), cf. 

When the test charge distribution is identical to that 
of the leading charge, t(r) = f(r), Eq. l(TH)l gives the 
discontinuity impedance for the given transverse charge 
distribution f(r) of the beam as 

Z(u) = -iZ - [e v {f; k)} 2 (a m + , (23) 
c p z 

where we use the effective polarizabilities a rn = Vv-r/2 
and a e = — x/2, and e„(/; k) is defined by (|19|l . 

The generalized longitudinal impedance is usually de- 
fined with the leading and test point charges having dif- 
ferent transverse displacement from the chamber axis, s 
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and t, respectively 0,0]- In that case, Eq. i(TB|) can be 
rewritten as 



Z(s,t;u) 
where 



-iZ — e v (s ; K)e v (t; re)(a m + -£] 
c [3 Z ' 



e„(r ; k) = - 



g 



*2 



(24) 



(25) 



The last expression follows from Eq. (|19|l because f g — 
e g {r ) for a point charge with transverse displacement 
r from the chamber axis. The impedance (|24|l in- 
cludes higher multipole longitudinal impedances. The 
usual monopole longitudinal impedance is obtained from 
Eq. (|24() when the charges are on axis, s — > and t — > 0: 



Z(u>) = -iZ o -el(0; re)(a m + ^|) . 

1 and 7 



In the ultrarelativistic limit, (3 - 
normalized field e„(0; re) becomes 



P :^-^fc 3 - 2 e9 (o)v,e£, 



(26) 
oo, the 

(27) 



cf . (|21|l , and Eq. 1(26|) coincides with the known result for 
the impedance of a small discontinuity [lj, |2| • 

For some particular cases of simple cross sections 5, 
one can obtain explicit expressions of the normalized field 
(|19|l and the longitudinal impedance 1)26(1 . The derivation 
is presented in Appendices for circular and rectangular 
cross sections. For a circular cross section of radius b, 
()26J) takes form 



Z(lo) = —iZq 



lu a m + (3' 



~'"' e [Io(Kb)}- 2 . 



(28) 



c 4ir 2 b 2 

which coincides, up to notations, with the result in |10| . 
Here Io(x) is the modified Bessel function of the first 
kind. 

For a rectangular cross section a x 6, assuming that 
the hole is located on the side wall at x — a, y — y^, the 
longitudinal impedance ((2611 is 



„, v .„w« m + (3 2 a e 

Z(uj) = -iZ Q -2 

c b £ 



E 

.P=0 



(-l)Psin[7r(2p+lW^] 
cosh(7rw 2p+ i/2) 



(29) 



where u m = a^Jm 2 /b 2 + re 2 /ir 2 . In the ultrarelativistic 
limit, Eq. (|29|) coincides with the result in 0, see Ap- 
pendix B for detail. 

For the case of an axisymmetric small obstacle on the 
wall of a circular beam pipe — like a small enlargement 
(cavity) or an iris — the longitudinal impedance turns 
out to be very similar to Eq. (|2*6*)l : 

Z{tu) = -iZ^el(0; K )2nb{a m + ^) (30) 

U3 a m + [3' 2 a e 2 
= - lZo ~c 2^6 [Io{Kb)] ' 



where the effective polarizabilities a m and a m are now 
defined per unit length of the circumference 2irb of the 
chamber cross section (circle) S. In transition from the 
first line in Eq. H3U|) to the second one, we used the 
known expression for ej,(0; re) in a circular cross section, 
cf. Eq. 1(28(1 . Essentially, for an axisymmetric discontinu- 
ity a m ,e in Eqs. 1(26(1 and 1(28(1 are replaced by 2 / Kba m ,e- 
In the ultrarelativistic limit, 8 — > 1, Eq. <(3U(I coincides 
with the previous results 0, 0] . It is worth mentioning 
that if an axisymmetric enlargement has area A of the 
longitudinal cross section, its magnetic polarizability is 
a m = A, while for an axisymmetric protrusion (iris) of 
the same cross-section area a m = —A. The electric po- 
larizability a e can be found by solving a 2D electrostatic 
problem, see examples in 0,0, El- It is positive for pro- 
trusions and negative for enlargements, so that in both 
cases a m and a e have opposite signs. 

The longitudinal impedances ((24(1 . I(26|) . and also 13U(I . 
depend on the beam velocity in two ways: via e v {s; re) 
and in the combination of polarizabilities (a m + a e //3 2 ). 
The first dependence enters via the parameter re6 = 
w6/(/?7c), as one can see from Eqs. P% |l -l|5U )l . For 

Kb <C 1 the factor e 2 (s ; re) is close to its ultrarelativistic 
limit, while at re6 > 1 it decreases exponentially to zero. 
We should emphasize that for j3 < 1 the monopole longi- 
tudinal impedance depends on the beam position in the 
chamber cross section, unlike its ultrarelativistic coun- 
terpart . For a circular cross section this dependence 
takes a particularly simple form as an additional factor 
of I 2 (ret) in ((28(1 and ((30(1 . where t is the beam transverse 
displacement from the chamber axis, see Eq. I(A3(I in the 
Appendix. 

In the combination a m + a e / 1 the electric contribu- 
tion is enhanced as the beam velocity decreases. This 
sum (more exactly, this difference, because a m and a e al- 
ways have opposite signs) can either vanish for some val- 
ues of f3, or become much larger than its ultrarelativistic 
limit a m + a e . It vanishes when < /3 = y/ — a e /a m < 1. 
We should note that such situation occurs only for dis- 
continuities like holes or chamber enlargements (small 
cavities), since for them a m > \a e \ |2|. For instance, a 
circular hole of radius ft in a thin wall has a m = 4/i 3 /3 
and a e = —2h 3 /3, so that y/—a e /a m = l/y/2. The 
impedance ((26(1 of the hole, which is inductive for rela- 
tivistic beams, changes its sign for (3 < l/v2 becoming 
a "negative inductance" . On the other hand, for pro- 
trusions and irises the impedance remains inductive for 
any beam velocity because of a e > \a m \, cf. |E 0- As 
a simple example, a semi-spherical protrusion (bump) of 
radius a on the wall has polarizabilities a e = 2ira 3 and 
u m = -7ra 3 0. 

For all small discontinuities the impedance vanishes at 
very slow beam velocities, when (3 — ► 0, since the fast 
decrease of the factor e 2 (s ; re) suppresses the growth due 
to a e /f3 2 . This behavior is illustrated in Fig. 1 for a hole 
and in Fig. 2 for a protrusion. The impedance magnitude 
can exceed the ultrarelativistic value many times. In fact, 
the ratio Z{f3)/Z{\) for tub/c = 0.1 in Fig. 1 reaches -83.3 
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at /3 = 0.062, and in Fig. 2 its maximum is 167.5, well 
outside the shown range. The extremes become even 
larger for lower frequencies. 




FIG. 1: The ratio of the longitudinal impedance 128^ to its 
relativistic value for a circular hole in a round pipe versus /3 = 
v/c for ub/c — 0.1,0.25,0.5,1 (solid, short-dashed, dashed, 
and long-dashed curves). 




FIG. 2: The same, for a small semi-spherical protrusion. 

Another interesting example of /3-dependence is for 
long elliptic slots parallel to the chamber axis. If the el- 
lipse semi- axes w, I satisfy w <C I <C b, the leading terms 
(cx w 2 l) of a m and a e for a thin wall cancel each other 
0,0 in the ultrarelativistic limit: 

7TU- 4 ., 41 . 

a m + a e « -r7-(ln 1) . 

61 w 

For (3 < 1 there is no such cancellation, so that 

1 3/3\ 



TTW 2 l TTW 4 1+P 2 41 



which for small (3 has an opposite sign and larger mag- 
nitude compared to its relativistic limit. 

The longitudinal impedance in Eqs. (|23[l - i|26|l is purely 
inductive. In the sense of a perturbative expansion in 
the small parameters |a m |/& 3 <C 1 and |a e |/6 3 <C 1, it is 



the first-order term. In Ref. pj higher-order terms were 
studied by taking into account radiative corrections to 
the fields near the discontinuity. The second-order term 
has both an imaginary and real part; the last one appears 
only at frequencies above the chamber cutoff. The real 
part of the impedance includes contributions from only 
a finite number of the eigenmodes propagating in the 
chamber at a given frequency, i.e. those with k g < lu/c 
or k' g < lo/c. The dependence of Re Z on frequency is 
complicated: it has sharp peaks near the cutoffs of all 
propagating eigenmodes of the chamber, but increases 
on average with the frequency increase. Well above the 
chamber cutoff, i.e. when u)b/c 3> 1 (still ujfi/c <C 1 to 
justify the Bethe approach), this averaged dependence 
can be derived in exactly the same way as in [jj] for ul- 
trarelativistic beams. The result is 



Re Z(u) = f^e£(0; n){a 2 n + a 2 



(31) 



where a ms = cf. Eqs. 0-©- Alternatively, the 

same answer can be obtained by calculating the energy 
radiated by the effective dipoles Q into a half-space. The 
physical reason for this coincidence is that at frequencies 
well above the cutoff they radiate into the waveguide the 
same energy as into an open half-space. 



B. Transverse Impedance 

Since we are interested in the dipole transverse 
impedance, it is convenient to consider the excitation by 
a dipole: two point charges with the opposite signs dis- 
placed from the chamber axis by s and — s, correspond- 
ingly, instead of Eq. QJ, moving at velocity (3c along 
the axis. Such a change would modify the expression 
(12411 for the longitudinal impedance replacing e u (s ; k) by 
e* p (0;/«) = s- S7e u (0;n), which is the second (dipole) 
term of its Taylor expansion, with e^O; re) being the first 
(monopole) one. After normalizing the excitation to the 
unit transverse beam displacement (dividing by s), we in- 
tegrate the synchronous harmonic of the transverse force 
acting on a test charge displaced from the axis by t, and 
then divide by the harmonic amplitude of the dipole mo- 
ment creating the deflecting force: 



ZAuS) = —— / dz exp ( — i— I x 



(32) 



x E±(t ,z;uj) + Zofiz x H±(t,z;u) . 

where the usual definition of the dipole transverse 
impedance assumes the limit of t — * s — ► 0, e.g. or 
|17|. The integration result includes contributions from 
both TM- and TE-modes excited by the effective dipoles 
in the beam pipe, unlike the expression for the longitu- 
dinal impedance (|24|l . where only TM-modes contribute. 

There is an alternative way to find the transverse 
impedance. Since we already have expression (I24|) for 
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the generalized longitudinal impedance Z(s, t;w), we can 
apply the Panofsky-Wenzel theorem, e.g., in [TsJ . Ac- 
cording to the theorem, the transverse impedance can be 
derived as 



Zj_(s,t;uj) = — V t Z(s,t;uj) 

UJS 



(33) 



where the longitudinal impedance is calculated with the 
dipole excitation as discussed above. If we follow this 
way, which is simpler in our case, the result will include 
only TM-mode contributions, or, in other words, the EFs 
of Q , since only they enter Eq. Q24f) for the longitudinal 
impedance. We used the relations between the trans- 
verse components of the beam fields, cf. Eq. JSJ, to prove 
that the two results obtained are equivalent. One should 
remark that definition (|32|l sometimes includes an extra 
factor 1//3, e.g. in [17|. For consistency with the rela- 
tion , here we use the transverse impedance definition 
without that extra factor, i.e., Eq. (|32l) . 

The final expression for the transverse dipole 
impedance of a discontinuity is 



Z±(v) = -iZ (3(a m + —)[s ■ d(K)}d(K), 
P 



(34) 



where s = s/s is a unit vector in the direction of the 
beam deflection from the chamber axis, and 



d{n) = Ve„(0;/c) = - 



Ve 3 (0)V^ 



9 



U2 

K g 



(35) 



The impedance dependence on the discontinuity shape 
is obviously the same as for the longitudinal impedance. 
The direction of the vector of the transverse impedance 
(|34|l gives the direction of the deflecting force acting on 
a displaced beam. As one can see from Eq. ()34[) . this 
direction is defined by vector d, Eq. 135|) , while the force 
magnitude varies depending on the relative direction of 
the beam displacement s with respect to d via the scalar 
product s ■ d. 

For a circular beam pipe, the transverse impedance 
(|34|l can be written, making use of Eq. (|A5|) of Appendix, 
simply as 



Zj_(uj) = -iZ /3 



7T 2 6 4 



id) 



2h(Kb) 



(s-h)h, (36) 



where h = b/b is a unit vector in the chamber cross sec- 
tion S directed from the axis to the hole (discontinuity). 
In this case, one can rewrite the dot product in a more 
conventional form, s ■ h — cos(</? s — ^Ph), where ip s is the 
azimuthal angle of the beam position in the cross-section 
plane, and iph is the azimuthal angle in the direction to 
the hole. Therefore, the deflecting force in the circular 
pipe is directed to (or opposite to) the discontinuity, and 
its magnitude depends on the angle between the hole and 
the beam transverse displacement, as was pointed out al- 
ready in the relativistic case 0. Equation (EH agrees, 
up to notations, with the result obtained in [10|. 



In a general case, it is sometimes convenient to rewrite 
the dipole transverse impedance (|34|l as 



Z±(lj) = -iZ Q (3(a m + -|) d 2 (K)dcos(ip 



9d) 



(37) 



In particular, this form is more convenient for the rect- 



angular chamber. Here d = yd x + dy, where x,y are 

the horizontal and vertical coordinates in the chamber 
cross section 5, and d x and d y are projections of c?, cf. 
Appendix B; tp s = (ft is again the azimuthal angle of the 
beam position in the cross-section plane; d = d/disa, unit 
vector in this plane in the direction of d, and tfd is the 
corresponding azimuthal angle. Obviously, d x — d cos ipd 
and d y — rising- It is seen from Eq. (|37|l that the an- 
gle ipd shows the direction of the transverse-impedance 
vector Z± and, therefore, of the beam-deflecting force. 
Moreover, the magnitude of Zj_ is maximal when the 
beam is deflected along this direction and vanishes when 
the beam offset is perpendicular to it. For a circular pipe, 
= <Ph, as was discussed above. For a general cross sec- 
tion, this is not the case even in the relativistic limit; see 
[2| for rectangular and elliptic chambers. As an illustra- 
tion, in Fig. 3 we plot the relation between tp^ and iph 
for a square chamber for a few different values of param- 
eter nb — ojb/((3jc). The case Kb = corresponds to the 
relativistic (or low-frequency) limit; as Kb increases, the 
difference between Lp<i and <ph becomes smaller. 




20 30 

<Ph, deg 



FIG. 3: Direction ipd of the transverse impedance l|37^ ver- 
sus hole position on the side wall of a square chamber for 
ujb/{n(5'yc) = 0, 3, 10 (short-dashed, dashed, and long-dashed 
curves). The solid line is for a circular pipe. 

In Fig. 4, the behavior of c?(k) in Eq. (|37|l is illustrated 
for a square cross section of the vacuum chamber. This 
function is given by Eq. (|B6|) in the Appendix; for plot- 
ting d(K) is multiplied by ab to make it dimensionless. For 
comparison, the same plot shows a similar dependence for 
the circular cross section, see Eq. Figure 4 demon- 

strates also the impedance dependence on the position of 
the discontinuity on the chamber wall: yh/b = 0.5 corre- 
sponds to the middle of the wall, 0.75 is one quarter from 
the corner, and 0.9 is close to the corner. Figures 3-4, as 
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1 2 3 4 5 6 

z=cjb/(nfiyc) 

FIG. 4: Normalized magnitude of function d in Eq. (I37H ver- 
sus nb/n for three hole positions h = yu/b = 0.5,0.75,0.9 
on the side wall of a square chamber (short-dashed, dashed, 
and long-dashed curves). The solid line shows for comparison 
Kb/2/h(K,b) from Eq. lESl . 



well as Figs. 1-2, were produced using Mathematica |20| . 
It takes less than ten terms in the series <|B6(1 and IIB3I) 
to get a very accurate answer for the sum, which is not 
surprising since the 10th term in i|B3|) is of the order of 
0(e~ 30 ) = O(10~ 13 ). However, even with hundreds of 
terms Mathematica gives the results almost as quickly as 
those involving modified Bessel functions for the circular 
pipe in Figs. 1-2. 

One should remind that the impedance magnitude de- 
creases even faster with the frequency increase for a 
fixed j3 than shown in Fig. 4, since Z± oc d 2 . The full 
/?-dependence of the transverse impedance includes, of 
course, the shape factor (3a m + a e //3; overall, it is sim- 
ilar to that of the longitudinal impedance illustrated in 
Figs. 1-2. 



IV. SPACE-CHARGE IMPEDANCE 

We will show in this Section that the same approach 
as above — based on field expansions in cross-section 
eigenfunctions — works equally well for calculating the 
space-charge impedance of a uniform waveguide with an 
arbitrary simply-connected cross section. While this sub- 
ject seems to be somewhat off topic for the present paper, 
the result comes as a by-product of the considerations 
in the previous Sections. The synchronous harmonic of 
the longitudinal beam field produced in a uniform beam 
pipe with the cross section S by the current @ is given 
by Eq. ©. Using the impedance definition l|16fl with the 
finite integration length L and assuming the same trans- 
verse charge distribution for the test charge as for the 
source, t(f ) = f(r), we obtain the space-charge longitu- 
dinal impedance per unit length of the chamber 



for a beam with the transverse charge distribution f(r ). 
The space-charge impedance (|3*S|) is expressed in terms of 
eigenvalues k g of the boundary problem (@J and the ex- 
pansion coefficients f g of the charge distribution in EFs 
0' fg = fsdrf( , r)eg(r). Every term of the series in 
Eq. (|38|l is positive, so that the space-charge impedance 
always remains a "negative inductance" . It is important 
to emphasize that the dependence on the beam charge 
distribution is an essential feature of the space-charge 
impedance. For example, one can not use a pencil beam 
(point charge) in l)38|l : the space-charge impedance di- 
verges as the beam transverse size vanishes, e.g. [Ig. On 
the contrary, the geometrical impedances discussed in the 
previous sections are generally independent of the beam 
properties, at least in the relativistic limit, and for that 
reason they are usually calculated in the simplest way, 
that is with an on-axis pencil beam. 

Let us compare Eq. (|38|) with the conventional form of 
the space-charge impedance in the long-wavelength ap- 
proximation, nb <C 1, 



Z sc {uj) 



LU Zq 1 

~c /? V 2^ 



9l 



(39) 



where g^ is called the longitudinal g-factor, e.g. |18, 21] . 
One can conclude that the sum in l|38[) is proportional to 
the g-factor: 



s(/;«) = E 



i 2 



2tt 



9l(k) , 



(40) 



which now depends on frequency and beam velocity, in 
addition to its familiar dependence on the beam charge 
distribution. The usual value of ^-factor is obtained as 
the limit at k — > 0, i.e. g^ = 5l(0). Using a term-by- 
term consideration of series (|40f) . it is easy to prove that 
g-factor is positive, gz,(«) > 0, and that frequency cor- 
rections reduce its value: (?l(k) < <?l(0). 

In the well-studied particular case of a circular cross 
section, we consider the beam in the form of a thin ax- 
isymmetric ring of radius a, f(r) = S(r — a)/(2vra). The 
summation is performed in Appendix A, and the result 



S{f;K) 



Io(na) 
2vr 



Ko(na) — Io(ko) 



(41) 



I (Kb) 

coincides with the familiar one, e.g. in ^l|- I n the long- 
wavelength limit k6 < 1, Eq. I|41|l gives the well-known 
g-factor of a hollow beam S — In (6/a)/(27r). 

For a rectangular chamber of width a and height b, to 
simplify calculations, we choose somewhat exotic beam 
charge distribution: a centered hollow beam with a sim- 
ilar rectangular profile £a, £6, where the beam-size scale 
< C < 1- Even in a simpler case of a square chamber, 
a = b, the space-charge factor (|40|) of the square hollow 
beam (a x (a is given by a rather long expression (see 
Appendix B): 



Z sc (lu) 



c /3 2 7" 



E 



ki 



(38) 



7T 3 C 2 



E 

p=0 



sin[7r(2p + l)C/2] 
2p+l 
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^ sinh[7TO 2p+1 (l-C)/2] x 

u 2p +i cosh(7ru 2p +i/2) 

'sin[7r(2p+l)C/2] 
x < ^4: . cosh(7rM 2 p+iC/2)- 



(42) 



2p+l 
cos[7r(2p + l)C/2] 

U2p+1 



sinh(7ru 2p+ iC/2) 



where u n = ^/n 2 + (nct/ir) 2 . We found analytically that 
in the relativistic (or long- wavelength) limit, na — ► 0, 
the leading term of Eq. |@I5J is In (1/C)/(2tt), exactly the 
same as in the case of a circular pipe above. The se- 
ries converges slower than the series in (|29[) or (|B6t , 
so up to 2500 terms were included in our Mathematica 
computations for this series. Still it took only a couple 
of minutes of CPU time on a PC to produce Figs. 5-6. 
They illustrate the behavior of S(f; k) as a function 
of the frequency and beam size. It is interesting that the 
sum S(f, 0) (relativistic or static limit) is very close to its 
In asymptotic even for the beam size as large as one half 
of the pipe size, £ = 0.5, cf. Fig. 5. The frequency de- 
pendence of the space-charge impedance (more precisely, 
of the g-factor) plotted in Fig. 6 is weaker than that for 
the impedances of small discontinuities in Fig. 4. 




FIG. 5: Function S I42|l for a square chamber versus the ratio 
£ of the beam size to the chamber size a at three different 
values of w — na/n = 0, 1, 10 (short-dashed, dashed, and 
long-dashed curves). The solid curve shows the analytical 
relativistic limit ln(l/£)/(27r). 



V. DISCUSSION 

The approach of Refs. 0, 0, where the impedances 
of small discontinuities were calculated in the ultrarela- 
tivistic case, was extended to beams with an arbitrary 
velocity and transverse charge distribution. The analyt- 
ical approach presented above provides a general picture 
of the coupling impedances for small discontinuities of the 
vacuum chamber with an arbitrary cross section in a wide 
frequency range, up to frequencies well above the cutoff. 
The upper limit on the frequency is imposed by the ap- 



CO 



1 

0.8 
0.6 
0.4 
0.2 



2 4 6 8 

W=6L>a/(7T/JyC) 

FIG. 6: The same versus w = na/n for three different beam 
sizes £a, C = 0.001,0.01,0.1 (short-dashed, dashed, and long- 
dashed curves). 



plicability of the Bethe theory: the wavelength must be 
large compared to the typical size of the discontinuity. 

We concentrated mostly on the leading (imaginary) 
part of the impedances created by discontinuities. In a 
general case, it was shown that the coupling impedances 
of small discontinuities depend on the beam velocity in 
two ways: first, through the combination of polariz- 
abilities (a m + a e //3 2 ), and second, via the parameter 
k = Lo/((3~fc) that enters into the factors e„(re) for the 
longitudinal impedance, see Eqs. (|2"51) - (|3U|) . and d{n) for 
the transverse one, Eqs. I|34 |) -l|57 |l . The first factor con- 
tains all the impedance dependence on the discontinuity 
shape, while the second one accounts for effects of the 
chamber cross section and the discontinuity location on 
the cross-section boundary. The second factor behaves 
in the same way for various cross sections: it decreases 
from its value at k = 0, which can be interpreted as ei- 
ther the ultrarelativistic or long- wavelength limit, as k 
increases. The decrease rate depends to some extent on 
the cross-section shape and the discontinuity position, 
but for Kb > 1 the decrease is fast (exponential), as evi- 
denced by Eqs. (gSJ, (J3UJ) and l|5fjjl. as well as by Fig. 4. 
The overall impedance dependence on the beam velocity 
for small discontinuities was discussed in Sec. III. The 
most interesting feature is that the impedance magni- 
tude for some /3 < 1 can exceed its relativistic value many 
times, see Figs. 1-2. 

The impedances for the circular and rectangular cross 
sections of the vacuum chamber were derived from the 
obtained formulas by substituting corresponding EFs. 
The expressions can be useful for calculating the beam 
coupling impedances in particle accelerators with non- 
ultrarelativistic beams. For the circular pipe our results 
agree with those obtained earlier for (3 < 1 in [lfl Hll | . 
In a similar way, the results for an elliptic chamber can 
be expressed in terms of Mathieu functions; it seems un- 
likely, however, that such a series would be convenient 
for calculations. For a complicated cross-section shape 
the relativistic- limit factors e^O) and d(0) can be found 
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numerically by solving a 2D electrostatic problem, see 
for an elliptical pipe. Applying the results obtained 
above, one can give then an impedance estimate for a 
given beam velocity in such a chamber. 

We also demonstrated that the same technique — us- 
ing the field expansion into a series of cross-section eigen- 
functions - works well for calculating the space-charge 
impedance of a uniform vacuum chamber with an arbi- 
trary cross section. A generalized expression g(n) for the 



longitudinal space-charge (/-factor is derived. It depends 
on the beam velocity and on frequency, in addition to 
the usual dependence on the transverse beam-charge dis- 
tribution. The generalized g-factor g(n) decreases mono- 
tonically as n increases: g(Q) > g(n) > 0. It is worthwhile 
to notice that the long-wavelength (or relativistic) limit- 
ing value g(0) can be calculated numerically even for very 
complicated vacuum chambers, e.g., those with screening 
wires and ceramic insertions 12111 . 
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APPENDIX A: CIRCULAR CHAMBER 

For a circular cross section of radius b the eigenval- 
ues k nm = [J. n rn/b, where /i nm is roth zero of the Bessel 
function J n (x), n — 0, 1, 2, . . ., and m = 1,2,.. .. The 



normalized EFs in circular co-ordinates (r, ip) are 

Jn(knm,r) J COS flip 



E 

nrn 



sin nip 



(Al) 



with N^ m = Trb 2 e n J 2 +1 (u nrn )/2, where e = 2 and e„ = 1 
for n ^ 0. For TE-modes, k' nm = fJ,' nm /b with J' n {n' nm ) = 
0, and 



i(r, if) 



Jn{k' nm r) ( cos nip 
,/N H 1 sinnip 



(A2) 



where N% n = wb 2 e n (l - n 2 / fi ,2 m )J 2 (u' nm )/2. 

The series l|25|) for the circular cross section can be 
easily summed 



e„(f; k) 



2 ^ cos [n(ip - iph)] x 



n=0 



E 



f-t'-nrn J n (/^nm^*/^) 

+ {Kb) 2 



J_ y> cos [n(ip - ip h )] I n {nr) 
TT&f^ e„ I n (nb) 



, (A3) 



where I n (x) are the modified Bessel functions, and n — 
u>/{(3"fc). For the last step in the above equation the 
summation was performed using formulas from [19 ] . The 
impedances of a hole in a circular pipe for [3 < 1 were de- 
rived in |9( in terms of sums similar to that in the second 
line of Eq. (|A3|) . The resulting sum l|A3() is essentially 
a multipolc expansion. In the limit of an on-axis beam, 
r — > 0, it becomes simply 



e„(0; k) = 



2nbI (Kb) 



(A4) 



The result simplifies even further in the relativistic limit. 
For 7 — + oo, e„(0; k) — > e„(0; 0) = e„ and becomes e v = 
l/(27r6), which also follows from the Gauss law. Then 
the inductive impedance 126fl takes an especially simple 
form derived earlier in Refs. |4],yj. 

For calculating the transverse impedance, we need to 
find the gradient of the field i|A3|l at the origin, cf. defi- 
nition J2HJ- Only the second term (n = 1) in the series 
131) gives a non- vanishing contribution when r — > 0: 



d(n) = Ve„(0;K) 



1 



nb 



nb 2 27i(k6) 



(A5) 
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where h = b/b is a unit vector in the chamber cross sec- 
tion S directed from the axis to the hole (discontinuity). 
In the limit of 7 — > 00, Eq. (|A5|) simplifies to d = h/(irb 2 ), 
in agreement with [j, O4] • 

For space-charge impedance calculations, we consider 
the beam charge distribution in the form of a thin ax- 
isymmetric ring of radius a, f(r) = S(r — a)/(2ira). Then 
the EF expansion coefficients are 



fn 



Jo(Ho m a/b) ( 1 



(A6) 



where 5 nm is the Kronecker symbol. The sum <|4U|) be- 
comes 



1 

S(f; K ) = - ^ 



JoifJ-Oma/b) 



W ^1 J l(V0m) (Mom + K 2 b 2 ) 



2ir 



Ko(na) — Io(na) 



K (nb) 

I {Kb) 



,(A7) 



where formulas |0] were applied. 



APPENDIX B: RECTANGULAR CHAMBER 

For a rectangular chamber of width a and height b the 
eigenvalues are k nm = n^/ n 2 /a 2 + m 2 /b 2 with n, m = 
1,2,..., and the normalized EFs are 



2 . irnx , irmy 
sin sin — ; — 



nb 



(Bl) 



with < x < a and < y < b. Let a hole be located 
in the side wall at x = a, y — yh- Assuming the beam 
displacement s from the chamber axis (a/2, 6/2), we can 
perform one summation in Eq. (|25|l to reduce the double 
sum into a fast-converging series 



where 



. p=° n, 
which coincides with the result 1] for the rectangular 

chamber. For some particular values of v, e.g. v = 1/2, 
the sum £(it,i>) can be expressed in terms of the com- 
plete elliptic integrals; in general, it is easy to calculate 
the series numerically because of its fast (exponential) 
convergence. Its behavior versus v for different values of 
the aspect ratio u was plotted in Ref. @|. 



(-l)Psin[7r(2p+ l)v] 
cosh[vr(2p+ l)u/2] 



(B5) 



The gradient of the field (|B2|I at the origin is required 
for calculating the transverse impedance. In the Carte- 
sian co-ordinates {x,y}, it is 



d{n) = {d Xl d y } — Ve„(0; k) = 

jy\ 1 N P u 2p+isin[7r(2p+ l)yh/b] 
ab \p~o sinh(7ra 2p+ i/2) 

2p sm(2Trpy h /b) \ 



cos 



h(7TM2p/2) J 



(B6) 



In the ultrarelativistic limit u m — > too/6, and gradient 
(|B6(I becomes frequency- independent: 



d(0) - 



- fvf lV , (2p+l)sin[7r(2p+ 1)^/6] 
b 2 \^ [ > smh[n(p + l/2)a/b} : 

2p sm(2irpy h /b) \ 
^ ' cosh{irpa/b) J ' 1 



e v (s; n) 



- 2^ sm[7rm(- + -^)] sin(— — ) x 

m=l 

sinh[7TM m (l/2 + s x /a)} 



sinh(7rw m ) 



(B2) 



For calculating the space-charge impedance of a cen- 
tered rectangular hollow beam £a x £6 in the beam pipe 
with rectangular cross section a x 6, the EF-expansion 
coefficients are 



where u m = ^/w?a 2 jb 2 + K 2 a 2 /-7r 2 . For an on-axis beam, 
s — > 0, it can be simplified to 

, n , l^ (~iysm[7r(2p+l)y h /b] 

e v (0;K) = - > j- . (B3) 

" p~o cosh[7ra 2p +i/2J 



In the relativistic limit, /3 — ► 1, 7 — *■ 00, Eq. (|B3I) becomes 
e,(0; k) -> e v (0; 0) = g„ = (~ ^) , (B4) 



fnm — 



4 sin(7rn/2) sin(7rm/2) 



2 C(a + 6)Va5 



(B8) 



a 7rnC 7rmC 

x [ — sin cos 

n 2 2 



b irnC irmC 

— cos cos 

to 2 2 



Substituting f nm into Eq. (|40|l and performing one sum- 
mation to reduce the resulting double sum to a series 
leads to Eq. (pE2"|) . 



